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ABSTRACT 


Expert sytems are? probably the most useful contribution of 
Artificial Intelligence. A number of expert systems are presently 
in use worldwide and many of them have come up with flying 
colors; but very little is done for developing expert systems in 
the field of numerical computations. 

Solving Partial Differential Equations (PDEs) numerically is 
indeed a complex process and needs a lot of expertise for 
choosing numerical & computational scheme, handling error control 
and stability etc. Moreover writing and debugging special purpose 
programs in procedural languages male the complete process 
cumbersome and trouble shooting. 

An expert system MA, is proposed and implemented, which can solve 
quasi-1 inear , second order PDEs. MA treats the problem in a 
generic way and addresses the physical problem directly. tt is 
developed in a relational environment, Prolog, which does the 
problem identification, preprocessing and postprocessing work. 
Since relational languages are inefficient for performing 
procedural jobs numerical computation is done in a procedural 
module, written in C. MA directly calls the procedural module 
’intelligently*, instead of generating procedural codes. MA is 
provided with a graphics module, which gives a graphical output 
of the numerical computation. Results are stored in a database 
chain and pointers to every set of data are stored in a B+ tree 
for efficient handling of the data. MA has a knowledge 
acquisition module, which allows the user to update the knowledge 
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base 


MA provides the user a very high level of interface and 


demands minimal knowledge about the problem in hand. An on 
help file is available for guiding users? moreover the user 
always ask why questions for clarifying his doubts. 

The implementation of the proposed system has been done 
PC/AT/XT using Turbo Prolog and Turbo C. The system operates 
MS DOS environment. 


line 

can 

on 

in 
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Chapter one - Int roduct ion and Review of Object Expertise 

1.1 Introduction 

With the increase in the computational ability of machines, the 
expectations from a user friendly scientific software is also 
increasing. The nature of the presently available software for 
computational problems has been rattier typical - a library of 
FORTRAN subroutines, which follow a srictly algorthmic path and 
obviously cannot process its ’'knowledge* to solve user's problem. 

Plethora of numerical schemes along with their stability 
and error propagation properties can often confuse the user. 
Moreover partial differential equat i ons ( PDEs ) are grouped 
into several complex hierar cliical status, each of which needs 
specific treatment to solve. A lot of expertise is needed to 
develop a lengthy code in procedural languages ( eg . Fo rt ran , Basic, 
Pascal etc.) for everyday applications in engineering design, 
fluid flow analysis, chemical engineering problems, heat transfer 
problems etc. Undoubtedly there is a serious need for encoding 
the PDE solving expertise. A system with a knowledge base 
containing this expertise will be able to select the most 

suitable problem solving strategy. 

A number of research groups started working with the above goal. 

In the past few years Problem Salving Environment systems for 

PDEs have been developed C1-43. They generally take the form of a 

procedural code generator along with a higher level user 

interface . 

MA is an expert system for solving PDF oriented problems. Fig.1 
depicts its overall logical structure. Present version of MA can 
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handle second order , linear or nonlinear PDEs. Instead of 
generating a code in procedural languages for numerical 
computation, MA ’intelligently’ calls procedural number crunching 
modules, whenever needed. MA’s conception of PDEs and numerical 
solution methodology is realized in terms of an object oriented 
representation which is implemented in logic programming 
env i ronment . 

1.2 Brief survey of the available Problem solving environments 
for partial differential equations 

The concept of automated problem solver evolved in late 70* s. 
Initially the research was directed towards the sole use of logic 
programming languages; but later, the inefficiency of these 
languages for procedural tasks was exposed. From early 80* s 
automatic. programming environments sorted coming up C1-S7J. Some 
of the systems start with continuous description of a problem and 
others work from the PDEs which are drawn from the conservation 
laws. Finally, some systems use symbolic manipulation to derive 
all or parts of the discretization scheme while others use 
artificial intelligence techniques to classify the type of the 
problem being solved and to fill program templates which 
implements the discretization scheme. Although a number of good 
automatic programming environments are available, the major 
drawback of most of them is the appl i cabi 1 i ty for a narrow class 
of problem; this is because these systems have knowledge about 
specific kinds of problem without having the deeper 

interpretation of the knowledge. It is indeed very difficult to 
go beyond these special purpose solvers using the conventional 
techniques of software engineering. Above all, most of them 
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suffer from serious problems about user interface. Complex 
fields, like numerical solution of PDEs , demand something more 
than a simple user friendly interface; user should be able to 
treat the problem and the knowledge in a symbolic way, which is m 
fact the paradigm of Artificial intelligence and Expert systems. 
Sytem should at first identify the problem physically and then 
only its mathematical interpretation in terms of PDEs should be 
made. System which start with discrete descriptions of the 
problems are likely to be more flexible but at the cost of the 

user interface, mentioned above; this also required enough user 
expertise because the user must deal with both the continuous 

pi r obi e m and i t s d i s c. r et e a pp>ro x i. mat i o n . 


i.3 Knowledge elicitation from the world of partial 
equations 


differential 


Human mind has been searching for the truth behind the 
mystic nature from the early days of its c i v i 1 i zat i on . It always 
tried to explore the nature by simulating it either theoretically 
or physically. The very birth of mathematics was endowed to the 
eternal search for a symbolic language, which would be able to 
describe the nature more precisely. Partial differential 
equation is a certain kind of mathematical expression, which is 
used for modeling several physical systems. 

Majority of the problems of physics and engineering fall 

natUraUy int ° °" e of the Physical categories: equlibrium 

problems , eigenvalue problems, and propagation problems. 
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Equilibrium problems are problems of steady state m which the 
equilibrium con F igurat i on 0 in a domain D is to be determined by 
solving the differential equation 

LC03=f — (1.3.1) 

within D, subject to certain boundary conditions 

B j r0D=g j — (1.3.2) 

on the boundary of D. Very often, but not always, the integration 
domain D is closed and bounded. In mathematical terminology such 
problems are known as boundary value problems. Typical physical 
examples include steady viscous flow, steady temperature 
distributions, equilibrium stresses in elastic structures, and 
steady voltage distributions. 

Eigenvalue problems may be thought of as extensions of the 
equilibrium problems wherein critical values of certain 
parameters are to be determined in addition to the corresponding 
steady-state conf igurations . Mathematically the problem is to 
find one or more constants (a), and the corresponding 
f un ct i ons ( 0 ) , such that the differential equation 

LC03=ocMC03 — (1.3.3) 

is satisfied within D and the boundary conditions 

B i C03=ocE i C0D — (1.3.4) 

hold on the boundary of D. Typical physical examples include 
buckling and stability of structures, resonance in electric 
circuits and acoustics, natural frequency problems in vibrations, 
and so on. 

.. arp initial value problems that have an 

Propagation problems are mniai 
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unsteady state or transient nature. Given the initial state, the 
subsequent behavior of a system can be predicted by solving the 
differential equation 

LC03«f (1.3.3) 

within the domain D when the initial state is prescribed as 

I j_C0!3=h i (1.3.6) 

and subject to prescribed conditions 

B i C03=g i (1.3.7) 

on the open boundaries. These problems are also known as initial 
boundary value problems. Typical physical examples include the 
propagation of pressure waves in a fluid, propagation of stresses 
and displacements in elastic systems, propagation of heat, and 
the development of self excited vibrations. 

Partial differential equations are classified according 

to their order, linearity, and boundary conditions. 

The order of a PDE is determined by the highest order partial 

derivative present in that equation. PDEs are categorized into 

linear, quasilinear, and nonlinear equations. Consider for 

example, the following second-order equation ' 

a(.)d£u + 2b( . ) d£u + c<.)d£u + d(.)==0 (1.3.S) 

dy s dxdy dx s 

If the coefficients are constants, or functions of the indepedent 
variables only C( . >2<x,y>3, then the equation is linear. If the 
coefficients are functions of the dependent variable and/or any 
of its derivatives of lower order than that of the differential 
equation C ( . )=( x ,y , u, du/dx , du/dy ) then the equation is 
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Neumann conditions Thu derivative oT the dependent variable 
is given as a constant or as a function of the independent 
variable. For example, 

dT/dx=G at x~1 and t >=0 

In a physical problem a neumann boundary condition can be 
specified if the gradient of the dependent variable is known, for 
example perfect thermal insulation in case of a heat conduction 
problem may give rise to the condition dT/dx-0 at the insulated 
boundary . 

Cauchy conditions A problem which combines both Diriclilet 
and Neumann conditions is said to have Cauchy conditions. 

On the basis of their initial and boundary conditions, PDEs may 
be furthar classified into initial value or boundary value 
problems. In the first case, at least one of the independent 
variables has an open region, in case of unsteady problems, time 
variable has the range 0< -t < =a and no condition is specified at 
t ~<x is an initial value problem. When the region is closed for 
all independent variables and conditions are specified at all 
boundaries, then the problem is of the boundary value type. 

1.3.2 Numerical solution 

The dynamics of physical systems that have more than one 
independent variable can be modelled by using PDEs. Only a few 
types of them render themselves to analytical solutions, for 
example some linear parabolic and elliptic PDEs. However, the 
majority of PDEs, especially the nonlinear ones require the 
application of numerical techniques for their solution. 
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Numerical solution of PDFs in itself is a wide field and a 


plethora of schemes exists For tackling them. No attempt would be 
made to describe all of them; only, the Finite difference scheme 
will be discussed briefly. 

1.3.3 Finite Difference Approximations 

The calculus of finite differences may be characterised 
as a “‘two-way street" which enables the user to take a 
differential equation and integrate it numerically by calculating 
the values of the function at a di s c rete ( f in i te ) number of points. 
In differential calculus, the definition of the derivative is 
g i ven as 

df < x ) /dx i !<0 =f ' (xo)--lt f (x)-f (xo) 1.3.10) 

x~>xo x - xo 

In the calculas of finite dif f er ences , the value of (x-xo) does not 
approach nero but remains a finite quantity. If we represent this 
quantity by h, 

h = x - xo 

then the derivative may be approx imated by 

f ’ ( xo f ( x)~f (xo) (1.3.11) 

h 

This concept of replacing the continuous operators by discrete 
ones can be extended for backward, forward and central difference 
approximat ions by using Taylor series as shown below: 

Development of Taylor series for u(x+5x,y) about (x,y> gives 

u(x+<Sx,y) = u(x,y) + &'.< du( v y > + ( <5x s /2 i ) d£u,( x y ) 

dx dx3 

+ ( (<$x> 3 /3! )d~M(x,y) +0<<5x 4 ) 
dx 3 
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which upon division by fix, results in the relation 

£Li±(x,y) ~ Fu(x-+<$x,y> •• u ( x , y ) U/tfx + 0<£x) --(1,3.12) 
dx 

This is known as the forward difference approbation of the first 
order differential operator. Similarly using the Taylor series 
for u < x-~<Sx , y ) about <x,y) backward and central difference schemes 
can bo obtained. 

Table 1 shows a listing of the Finite difference approximations 
of partial derivatives using central differences. 

1.3.4 Stability and convergence 

The presence of round-off error or any other 
computational error may lead to numerical instability. 
Consider the simple ordinary differential equation 

dy/dx = y - x — (1.3.13) 
whose general solution is 

y =- Ae :< + (x+1 ) --(1.3.14) 

where A is the integration constant. In an approbate solution 
the eponential term is likely to be introduced as a resul t of 
round-off errors. Any numerical scheme which allows the growth of 
error, eventually "swamping ' the true solution, is unstable. 
These numerical phenomena must be avoided by restrictive action, 
such as limiting the interval size or adopting an alternative 
method . 

Let us talk about stability in a more mathematical way. Let 
U(x,t) be the solution of a given difference appro bat ion , 
solvable step by step in the t direction. The effect of a 
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mistake or, more likely, round off error in machine computation 
may replace U(x 0 ,t 0 ) by U(x 0 ,t 0 > + e at the grid point (x 0 ,t 0 ). If 
the solution procedure is continued with the value U(x o ,t 0 > + e 
without new errors being introduced, and if at subsequent points 

•Jf* *$$• 

the value U (x,t)is obtained, then we denote by U (x,t)-U(x,t) 
the departure of the solution resulting from the error e at 
( x o , 1 0 ) . When errors are introduced at more than one point, 
cumulative departure's result which are not additive ecept in 
linear problems. If £ is the maximum absolute error — and h 
the interval size, then a procedure is said to be pointwise 
stable if the cumulative departure tends to zero as 5 — ..= 0 and 
does not increase faster than some power of h ^ as h — >0. 
when the corresponding continuous solution remains bounded a 
finite difference process within the semi-infinite strip CKx<1, 
t>G is called stepwise unstable if for a fied network and fied 
homogeneous boundary conditions there exists initial disturbances 
for which the finite difference solutions become unbounded as 
j — >« . 

A second fundamental concept, that of convergence, is often 
related to stability. To introduce this idea we utilie the PDE 
L(u)=0 in a region D, u=g on T 

where T is the boundary of D. Let us suppose that there is 
only once, say h, and write the finite difference problem as 
Lf,<U)=G in D, U=gj- t on T. 

We say the finite difference scheme converges if U (P ) 
converges to the solution u(P) , with the same boundary values, as 
h — >0. 
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1.3.5 Integration in time domain 

Integration of the dependent variable can be approximated by 
a weighted backward difference along the time axis. Depending 
upon the weighting factor several schmes are proposed. 

For example, a second-order partial derivative can be expressed 
as a weighted average of the central differences at points 
(i,j+ 1 ) and (i,j), where 3 denotes the spatial direction and j 
denotes time direction s 

d s u/dx s — ©<S s u j ? j + -j + ( 1 ~ 0 ) <9 s u 3 ^ j 

= ©IK '-H-H ,j-M " 2'J- j. , _H-i + u i- 1 , j + 1 >/h2'.1 + 

+ (1-0)C<u i + 1/ j - 2 u iyj + u :i i , j ) / h D 

- - (1.3.15) 

where 0 is in the range Q<0<.1. Similarly other differential 
operators can also be integrated with a time weighting factor. 

Explicit Scheme s 

When 0=0, equation 1,1 reduces to its second part, 
which involves property values at the previous time step. This is 
known as explicit scheme. The stability of this scheme is 
conditional . 

Crank-Nicolson Scheme ' 

When 0=1/3, equation 1.1 yields the Crank-IMi col son 
formula for second order differential operator. This is also 
conditionally stable. 

Implicit Scheme - 

Finally, when 0 = 1 , equation 1.1 becomes the implicit 
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formulation. Implicit formulas are un condi t ional ly stable. 

It can Le generalized that most explicit finite difference 
approximations are conditionally stable whereas most implicit 
approximations are unconditionally stable. The explicit methods, 
however, are computationally easier to solve than the implicit 
techniques . 


1 .4 Present work 

Several limitations of the present problem solving environment 
(PSE) for partial differential equations, as mentioned in section 
1.2 have motivated many researchers M-63 to work for systems 
with greater flexibility for handling equations along with a 
relational handling of the vast knowledge base. 

We propose a system, lift which has several unique features for- 
solving PDEs and equipped with some powerful AI techniques. MA 
eploits the properties of both relational and procedural 
languages. It has a robust PDE solver along with a higher level 
user interface. (, F i g . 1 . 3**4 shows a schematic diagram of the 
system. The main goal of the expert systems is to treat the 
problem in a natural symbolic way. MA aspires this by handling 
the PDE solving problem as a physical problem and is not limited 
to the computational procedures. 

MA consists of three main modules s 

1) Preprocessor 

2) Solver 

3) Postprocessor 

The logical organization of MA is shown in Fig.1. 
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Preprocessor module identifies the problem with the help of a 
physical problem knowledge base and receives required parameters 
for defining it uniquely. It also points out errors in typing. 
User defined parameters are checked for stability iriterias, 
specified in the rule base. 

Solver is a procedural module which does the numerical 
computation for solving the set of algebraic equations, generated 
from the discretized PDE, which is written in C. This module 
directly does the computation work by calling suitable number 
crunching procedures, instead of generating a code in procedural 
language and then making it an executable file, which is a 
conventional. feature for most of the automatic PDE solver. This 
unique feature makes MA more competent for microcomputer usage. 
Postprocessor module handles the solution kept in database chain. 
In case of time dependent problems times are kept in a B+ tree 
chain along with a pointer to the corresponding solution values. 
This makes the search for the solution at a specific time from a 
large set of solutions easier and fast. Adequate graphics 
facilities are available in the Postprocessor for plotting 
different kinds of graphs. 

MA treats the problem in a generic way and is not 
restricted by any ad hoc implementation strategy. This will be 
described in chapter 3 in details. 
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Chapter Two ~ Logic Programming Prolog & Expert Systems 

2.1 Introduction 

Logic programming is an approach to computer science in 
winch the horn clause form of the first order predicate logic is 
used as a high level programming language. The study of symbolic 
logic goes back to the work of Aristotle in the fourth century 
B.C. First order predicate logic is a branch of symbolic logic 
that has largely evolved in the twentieth century. Logic 
programming is based on a subset of first order predicate logic, 
but is equally broad in scope. Logic programming allows a 
programmer to describe a situation with formulas of predicate 
logic and then to use a mechanical problem solver to make 
inferences from the formulas. 

2.2 History of symbolic logic 

What i.s now known as "traditional logic" began at the time 
of Aristotle over 22 centuries ago. Aristotle attempted to codify 
into a scientific system the way that knowledge can be most 
effectively pursued through rational debate. 

Aristotle analysed the form of a statement into the following 
e 1 ements s 

Quantifier Subject Copula Predicate 

so that, for instance, in the following statements 

Some IITK students are intelligent 

"IITK students" is the subject, "Some" quantifies the subject, 
"are"is the copula and "intelligent" is the predicate. 



Traditional logic recognises four possible relationship between 
two classes. They ares 

All A is B. 

No A is B. 

Some A is B. 

Some A is not B . 

where A stands for a subject class of items and B stands for a 
predicate class of items. Hellenistic Greece also talked about 
syllogism. The syllogism is a set of rules which govern the 
conclusion from a set of premises. But Aristotelian logic was 
incomplete and inadequate to describe real world situations. It 
was the English mathematician Deliorgan and Boole who criticised 
Aristotelian logic first in the middle of the nineteenth century. 
Although traditional logic is a logic of classes, it contains no 
notion of complement. Moreover its complete dependence on natural 
language has led to many absurd difficulties. 

Boole and DeMorgan associated logic with mathematics rather than 
philosophy. They developed the concept of set theoretic operators 
for logic. That was the birth of modern symbolic logic which led 
to Propositional and Predicate logic. 

2.3 Propositional and Predicate logic 

2.3.1 Atomic 

Formulas 

A proposition is any statement that 
can be assigned a truth value (either true or false). In the 
conventional lexicon of propositional logic, atomic propositions 
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are defined 


P r °P»^t>on e that cannot bo broker, down into 
components. In Propositional logic, tl.e basic truth-valued object 
ii> an atomic proposition that cannot bo subdivided into 
components in any way. The basic truth valued object in predicate 
logic is an atomic formula. An atomic formula is composed of a 
predicate symbol, and terms that act as arguments to the 
predicate symbol. In general, the predicate symbol Is the name of 
a relationship that holds between the arguments. As an atomic 
formula is itself constructed of other objects, it posseses much 
greater expressive power than an atomic proposition in 
propositional logic. 

An atomic formula is written as a predicate symbol. followed fay 
some number of arguments inside parentheses. Each argument is a 
term. The general form of an atomic formula is s 
P(t1,t2 tn) 

where P is a predicate symbol and t1,t2.«...tn are terms. 

Terms 

A term can be one of the three things : 

1) a constant; 

2) a variable; or 

3) the application of a function. 

The application of a function is written as a function symbol 
followed by at list of arguments inside parentheses . Each argument 
is itself a term. The general form of the application of a 
function is ' 

f (t1,t2, tn) 

where f is a function symbol and tl , t2. . « . tn are terms. 
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2.3.2 Well formed 

formulas 

The result of combining atomic formula together with 
logical connectives is a well-formed formula <wff). A well formed 
formula can be defined as follows : 

1) an atomic formula is a wff? 

2) if A and E-i are wffs then so are 

Wff Read as 

'''A not A 


A & B 


A and B 


A v B 


A o r B 


A — > B 


A implies B 


A < : B 


A if and only if B 


for some 


Vx A 


for all 


Quant i f j. e rs 


and V are known as existential quantifier. 


As in propositional logic, the truth value of a wff 
depends on the interpretation .An inter pretation of a wff can only 
be made with respect to particular domain of inter pretat ion , 
which is the set of all possible values of terms that occur in 
the wff. 

2.3.3 Inferences in 

predicate calculus 

Inferences in the propositional calculus result from 
the internal structure of sentences that are comprised of 
component sentences. However, the inference in predicate calculus 
needs an analysis involving a finer breakdown of structure than 
one stopping at atomic sentences. Let us consider an example for 
a clear understanding of the subject-predicate relationship. 
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var i. bl e instantiation and existential quantifier. 

All kids are nice. 

A'l is a kid. <1> 

Hence, A1 is nice. 

Let kids(y) indicate the fact that the object y is a kid, and let 
niceCx) indicate the property that x is nice. Both kids< ) and 
nic.e( ) are called predicate of arity 1, since each involves one 
argument . 

using the notions thus far introduced, the inference (1) can be 
symbolized by 

V x C k i ds ( x ) — > niceCx) 3 
kids ( A1 > 

Hence, nicc(AI) 

Grammatically, the instantiated predicate nice(AI) represents the 
declarative statement "A1 is a kid." and is either true or false, 
verity is based on the distribution of facts. But what about the 
expression kids(x)? Factually, the expression is empty and is 
not a proposition. The set of objects that the variable in an 
open sentence might represent is called universe of discourse. 

2.4 Prolog - PROgramming in LOGic 

After discovering the computability of logical 
consequence in an abstract sense, it is a natural step to want to 
mechanize the proof process on a computer. Although full 
predicate logic is a very expressive language, implementation of 
proof procedures of predicate calculus in a general way may cause 
combinato r ial explosion. In 1965 Robinson came up with 
resolution, an inference rule appropriate for machine inference. 
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Loveland, Kowalski, and Kuehner -further refined this wnth the 


techniques of model elimination and the selection function. In 
early seventies, Colmurauer and Roussel brought this work 
together into a language that they called PROLOG (for PRGgramming 
in LOGi c ) . 

Following parts of this section describes resolution 
technique, which is the logical basis of the language, and 

several other features. 

2.4.2 Resolution 

Resolution worts as follows! 

two clauses can be resolved with one another if one of them 

contains a positive literal and the other contains a 
cor responding negative literal with the same predicate symbol and 
the same number of arguments, and if the arguments of both 

literal can be unified (matched) with one another. The clausal 
form of predicate logic is a way of writing formulas that uses 
only the connectives &, v, and <v . A literal is either a positive 
or negative formula. Consider a theory composed of the following 
two clauses: 

P (a) v *'Q(b, c) (1) 

Gt( b, c ) v '’“R ( b , c ) ( 2 ) 

Since clause (1) contains the negative literal '“G(b, c), and 

clause (2) contains a corresponding positive literal Q(b,c), and 
since the arguments of the two literals can be unified, then 
clause (1) can be resolved with clause (2). The resolvent is 
shown below • 

P (a) v ~R ( b , c) (3) 
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Now, all the three clauses can be used in future resolution. 

The basic problem is to prove whether a clause is or is not 
a consequence of a theory. To be able to automate the discovery 
process, we would like to find an efficient rule that can detect 
the inconsistency of a set of clauses. 

There are more than one problem solving strategy that can be 
pursued using the resolution rule. Prolog incorporates a top-down 
(or backward ) strategy. Tit is strategy aims to detect whether a 
single clause, C, is a consequence of an existing set of clauses, 
T.The set of clauses, T is assumed to be consistent. The 
algorithm works as follows! 

To begin with the negation of the clause is being tested, ~C , and 
is added to the existing set of clauses, to form a new set of 
clauses, T’.If the algorithm can derive the empty clause from T* 
it is inconsistent because of the presence of ' V C , and C must 
the r fore be a consequence of T. 

2.4.3 Prolog syntax 

The Prolog language is a combination of powerful ideas, 
including : 

1) the use of Horn clauses to represent knowledge; 

2) descriptive style of programming; 

3) both declarative and procedural semantics. 



In <-.i Horst clause 1 , one conclusion is followed by zero or more 

conditions, which is written in Prolog as follows : 

conclusion :~ 

condition 1 , 

conditions, 

conditionN. 

is read as "if", and , is read as "and"; so the whole clause 
can be read s 

The conclusion is true if 

conditionl and conditions and. . . conditionN 
are all true. 


2.4.4 Semantics of Prolog 

The semantics (meaning) of a formula of a symbolic 
logic. refers to its truth value. The semantics of a constant 
symbol inside of a formula refers to its value with respect to a 
domain of interpretation. In computing, on the other hand, the 
semantics of a programming language construct usually refers to 
the behavior of the computer when the construct is evaluated. 
Because Prolog is both logic and a programming language, both 
notions of semantics are applicable. 

There are three semantics models to explain the meaning 
of a Prolog program, which were first explained by Kowalski. They 
are : the declarative model, the procedural model, and the 
abstract machine model. 

The declarative semantic model of Prolog specifies the truth 
value of the relations. The word declarative is used because a 
Prolog clause declares that a relation holds between its 
arguments if all of the conditions of the clause are met. 


28 



For instance, the following clause : 


executive^ N«jh#e , Salary ) : - 

employee (Name, Salary) , Salary > 3000. 

can bo read as s 

Anyone (Name) is an executive if 
ho or she is an employee 
with a salary greater than Rs 3000. 

Read according to the declarative model, Prolog clause are 

formulas of first order predicate logic written in Horn clause 

form. 

According to the procedural model, the conditions of a clause 
specify a process to establish a truth value of the conclusion of 
the clause. A set of clauses with the same predicate name and the 
same number of arguments as a procedure is understood to be a 
call to that procedure. Read procedural ly, the meaning of the 
above clause is s 

One way to find an executive is : 
first, find an employee, 
then second, verify that the salary 
of the employee is greater than Rs.3000. 

A Prolog query has a computational meaning in the sense that they 
trigger a certain behavior on the part of the Prolog interpreter. 
The interpreter applies a problem solving strategy to evaluate a 
query against a set of clauses? its problem solving strategy can 
be characterized computationally by an abstract machine. The 
abstract machine model specifies the meaning of a query and set 
of clauses in terms of abstract machine actions. 
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2.5 Knowledge representation 

One of the main goals of the AI research community to js 
develop system which can simulate the real world knowledge in a 
machine. The art of encoding 'knowledge’ is technically known as 
knowledge representation. Wo have already discussed the 
representational formalism of first order logic. Now we shall 
discuss several means of representing knowledge briefly. 

2.5.1 Facts and Rules 

As we mentioned earlier, the a Prolog program consists of a set 
of clauses, which is either a fact or a rule. A fact is an 
assertion that a particular relation holds; it is written as a 
name followed by the arguments inside a parenthesis. The 

following fact expresses the idea that "ma solves pde" s 

solves (ma, pde ) . 

A rule is a fact whose truth value depends on the truth values of 
other facts. A rule can be represented in a Horn clause with each 
clause using assertions as well as the rules in the knowledge 
base. An example in Prolog is as follows 5 

X Y , Z. If Y and 7 , then X. 

X Y ; Z. If Y or Z , then X. 

2.5.2 Semantic network 

In a semantic network, entities and classes of entities are 
identified with nodes, and relations between entities are 
identified with arcs joining nodes. An arc connected to a single 
node establishes a property of that node. Several types of 
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attributes are often employed. For instance the class relation 
:is-a and the possesation relation has -a. In O-A-V 
(object -attribute --value) notation examples ofthis relations may 
loo k 1 ike » 


j. s__a ( pde , eq uat ion). 

is_a(ode, equation) . 

has__a( pde , two__inde pendent_var iables ) . 

An essential concept of the semantic network formalism is that of 
hierarchy, which makes it particularly adept at representing 
taxonomies of knowledge. Each level of the taxonomy is 
represented by a node connected by is„a arcs to higher and lower 
levels . 

2.5.3 Frames 

Minsky L13II originated the concept of a frame as a way to represent 
knowledge about situations. Each frame contains slots that 
identify the type of situation or specify the parameters of a 
particular situation. Moreover, other forms of knowledge, such as 
procedural information, can be used in slots. Inheritance is one 
of the important features of frame formalism. It is possible to 
specify that if a slot in one frame is not specifically filled, 
the frame will inherit a default value of that slot from a 
superior frame. Inheritance of slot values is enabled between two 
frames by the presence of an a_kind_of slot in one frame that is 
filled by the name of the other frame. The a_kind_of slot is 
similar to the is a arc in a semantic network. 
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E.6 Expert systems 


Expert systems are a new class of application programs, intended 
to make the knowledge of an expert in some special field readily 
available to "lay" users via a computer driven, interactive, 
dialog oriented system. The fundamental difference with respect 
to traditional programming is that the processing rules are 
stored in a knowledge base along with the data itself. Tho core 
of an expert system is an inference engine. It applies a 
particular strategy to draw conclusion from the knowledge stored, 
thereby producing new knowledge. Suitable heuristics must provide 
for targeted knowledge processing, so that a meaningful response 
to a user query can be produced within a reasonable amount of 
time. 

Fig. 2 shows a schematic diagram of a general expert system. 
Modern expert systems comprises two other important features 5 

1. knowledge acquisition module 

2. an explanatory interface 

The knowledge acquisition module can be seen as a kind of back- 
end to the knowledge base; and the explanatory interface as a 
front-end to the inference engine. But perhaps the most fruitful 
way of looking at this diagram is as a pointer to the four 
central unsolved problems in expert systems methodology . 

These are: 

1. the problem of knowledge representation 

2. the problem of approximate reasoning 

3. the problem of knowledge acquisition 

4. the problem of human-machine interfacing. 
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It may s e e m paradoxical to »ay that the state of art i.n expert 
systems is a collection of unsolved problems, when hundreds of 
practical expert systems are in worldwide, but it is true. 
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Chapter Three - System Analysis 


3.1 Introduction 

For the numerical solution of PDFs there is the never 
ending discussion of special purpose program versus general 
"black box" solver. An ideal Black: bo;; solver should have the 
following properties : Full flexibility of 

1) PDF. operator ( all types of linear *% nonlinear systems of 
PDFs) , 

2) boundary condition operator ( all types of nonlinear BCs), 

3) solution operator ( variable quality, error estimate, 
se 1 f adaptat ion ) , 

4) geometry ( arbitrary complicated 2-D or 3-D domain ). 

As we discussed in chapter 1.2 all available solvers are 
compromises. Our black box solver MA has the following properties 

1) solves quasilinear systems of PDFs, 

2) admits arbitrary nonlinear BCs, 

3) restricted to a rectangular domain (or a domain that can 
be transformed analytically to a rectangular domain ). 

Full geometric flexibility can only be obtained by FEM on 
unstructured grids, but then we have a data structure that 
ne ccess i tates much indirect addressing and thus reduces the 
efficiency of vecto r i zat i on . Above all we have for the FEM, 
problems with the requirements 1 to 3 mentioned above. The FEM 
thus needs always an individual adaptation to a distinct problem. 
This was the main reason for choosing FDM. 
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3.2 DOMAIN OF EXPERTISE 


Expert system MA was initially designed to serve as the 
computational front end of a knowledge based system for heat 
ex c hanger design. Even then , MA is flexible enough to treat its 
problem domain oP solving PDEs numerically in a generic way to 
other applications. 

At present it cam handle second order PDEs. This covers a 

wide range of engineering and scientific appl i cations .These 

equations can be expressed in the farm - 

A*d ? u/dx ? + B*d»u/dxdy -i- C*d*u/dys + D*du/dx + E*du/dy 

•v- F =•'• 0 -(3.2.1) 

where A , B , C , D , E »F and const are either constants or functions of 
any combination of x,y,u; u is the dependent variable and x,y are 
independent variables. Mon-linear terms are linearized by 
Jacobi's method using Taylor series. 

MA must be supplied with the required boundary conditions and 

possibly initial conditions in order to obtain a unique solution. 

It can handle boundary conditions in the form of Dirichlet » 

Neumann , Periodic or mixed conditions. After recognizing the 

problem with the help of its knowledge base it chooses an 

appropriate numerical algorithm and then goes for the solution. 

Knowledge base of MA contains facts and rules which it uses to 

achieve goals in the process of solving PDEs. It can always be 

updated and an experienced user can control the algorithmic path. 

Knowledge base rules contain the following : 

13 Classification of PDEs 
23 Convergence properties 

33 Applicability of di s c r et i zat i on schemes 
43 Choice of computing scheme 
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Some of the fact:; are stored in frames . Some of them are s 

•11 PDEs and their corresponding problem domains 
HI Discretisation scheme* 

31 Grid information 

41 Information about camputat s onal schemes 

3.3 KNOWLEDGE REPRESENTATION 

In order to reason about the algorithm chosen to solve the PDE » 
MA must have a knowledge base containing the expertise on problem 
domain . MA has an object oriented knowledge base about the world 
of PDFs , stored m rules and frames . We shall discuss the 
knowledge representation in detail here. 

Prolog implements a subset of first-order logic , the Horn 
clauses . A Prolog clause has both procedural and declarative 
interpretation . Constants , variables and structures can be used 
as arguments to predicates . More than one Horn clause may be 
needed to define a predicate and a corresponding AND/OR graph is 
thus created for each predicate. 

3.3.1 Rules 

The main clause of MA which controls the process of solving PDEs 
right from recognising it can be stated as 5 

start_process if 
get_problem(Equat ion , Type), 
get_gr id_parameters (Type , Param_l ist ) , 
disc:retise( Equation ,Param_l ist , Dis c r ete_l ist ) , 
sol ve (Discrete_l ist ) . 

start_process is a dummy predicate which fires a chain of actions 
, finally solving the PDE ; get„problem receives the problem from 
the user and then instantiates variables - Equation ( a list of 
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strings ) and Type (a structure containing the i dentif iet s oT the 
equation). These inherited terms are passed in 
g p t„_y r i d„pa r ame t e r s , d i s c f e t i z e an d s o 1 ve . 

3-3-2 Frames 

Frames are used in MA for storing straightf or wa r d informations. 
A frame is composed of slots for variables and their values. Here 
is a frame for storing discretisation schemes - 


D i s c r <? t i 7 a 1 1 o n s cl -i e m e 


n chert* el ! 

A 

( 

t 

schemes ! 

B 

t 

1 . .. * 

« M « 

1 1 

» t 


The corresponding representation in Prolog in given below s 

frame ( di scrot i sat i on__s theme , s charnel , a) . 
frame ( di s c r et i sat i on__s theme , schemes, h ) . 

Here schemel , schemes. . represent the names for the slots such as 
implicit , f ul 1 y_expl i cit etc. and a,b.. are structures containing 
cor res ponding approx imated< di s c rete ) terms. 

3.3.3 Object oriented representation 

Object oriented hierarchial represention has many remarkable 
advantages in contrast to the single stage classifier systems. MA 
classifies the problem from three aspects - mathematical » 
discretisation and computational . 

Fig. 3 depicts the hierarchical relationship among the objects 
based on mathematical taxonomy. A sequential search along the 
three main branches ~ equation , spatial domain and time domain 
finally identifies the object . Similar hierarchical relationship 
exists for the discretisation and computational schemes. 
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3.4 USER INTERFACE 


User interface is probably the most important factor for the 
survival of a software . liA has an interface designed to interact 
at a higher level through pull-downs 7 po p- up-men ues , editors 
and at the same time to be transparent enough to have a ’feeling* 
about the ongoing internal process . 

Fig. 4 depicts the main menu of MA. Two options are available 
for entering the problem — 

Feed equation User can directly feed the PDE to be solved 

through an editor. 

State problem s- Instead of explicitly mentioning the equation to 
be solved the continuous physics problem can also be stated , 
which in turn gathers necessary informations from the knowledge 

base . 

Default denotations of the dependent and independent variables 
can be changed by using option menu . 

Any error in typing the equation may cause unwanted errors in the 
computation. In order to prevent this , MA is provided with a 
smart parser which can point out the mistake. At any stage 
user can interrupt the system and enquire about the process. 

Plethora of discretization and computational schemes often 
confuse an -inexperienced user. MA always offer the best scheme to 
solve the PDE , but the user can also get rid of this. 

Graphics is an essential feature of MA. Its graphics modules has 
several features. 

Graphical output * 

Graphical representation of the result is often more 



comprehend ve than a simple set of numbers. MA is, provided with a 
graph plotting module which can also be called sc? paratel y . Fig. 3 
shows some graphical output of the system. Usual features for 
providing legends and lebels are available. At present the 
complete graphics module is written in Prolog. Handling and 
processing of solutions are unique and designed carefully for 
accessing them efficiently. All the solutions are kept in the 
form of external database, alongwith an index kept sorted in B+ 
trees with a pointer to the corresponding data. 

3.5 Controlling inference =- 

Backward chaining is one among the primary means of performing 
inference. In this scheme verity of a syllogism is checked by 
making certain that in each of the instances in which the 
conclusion is false at least one of the premises is also false. 
This is particularly suitable for our purpose , as the number of 
subgoals at every level of the hierarchical tree is limited. This 
process continues until either the goal is proven or failure 
occurs r in which case system backtracks to find alternate 
solutions. Object oriented hierarchy mainly controls the search 
path. All the properties of the superclasses are inherited. After 
defining the controlling path by traversing along the 
hierarchical tree, delicate control is performed by metalagical 
control. A sophisticated scheme for controlling the flow path at 
run time is implemented , which makes use of dynamic cut , 
assertion and retraction of informations about the problem such 
as type of equation , grid parameters etc. 
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3.6 Equation approximation and computation 


MA determines a] 1 the char acte r i st i. cs of an equation shown in 
Fig.1 by means of passing the equation recursively down through 
the three main branches of the hierarchy. Fig. 6 shows -some of the 
feed-backs to the user and the discretisation menu. liA solves 
PDFs by finite difference methods. For transforming the 
continuous domain PDF into a set of algebraic equations in a 
discrete domain , differential operators must be approximated 
MA offer'> several schemes for discretizing the equation , eg. 
Cranck -Ni colson , Fully explicit etc. for time integration and 
use central diff rente approximation for the spatial derivative 
terms. MA always checks for the convergency criteria of the 
chosen scheme against the grid parameters . It warns the user 
about the possible trouble and offers a chance to modify the 
parameters. If the user defines his own scheme by updating the 
knowledge base, he can always define the stability criteria for 
choosing appropriate grid parameters. 

Logic: programming environments are not that efficient to serve 

procedural jobs. MA incorporates a set of procedural functions 
written m C. Linear equations can be handled in a straight 
forward way; but for non-linear terms C-rnodules in turn call a 
symbolic computational module to simplify the terms. 
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Chapter Four Implementation 

4.1 Introduction 

So far we have discussed the basic knowledge of numerical 
solution of PDFs , concepts of Logic programming and a conceptual 
analysis of the proposed system, MA. In this section we shall 
talk about the practical implementation techniques and some 
special features of MA. 

4.P Why Prolog ? 

We have already seen some of the powerful features of 
Prolog in chapter II. Although Prolog is certainly not a very 
good language for procedural tasks such as solving simultaneous 
algebraic equations, matrix operations etc., it can be 
effectively used to develop a higher level interface for a 'lay' 
user. In fact it is the preprocessing part which needs an 
relational environment for identification and initial processing 
of the problem. Prolog is one among the best suited languages far 
this kind of job. We exploited the inherent backward chaining 
feature, backtracking and pattern matching to develop MA. 

Major problem with procedural languages is that the meaning of a 
construct is confined to itself and it is very difficult to 
develop a relational formalism among them, which is a must for 
identifying the problem in a hierarchical search space. 

4.3 Why Turbo Prolog ? 

Turbo Prolog is a comparatively younger version of 
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p i o . 1 . o g « 


i-.vcn 


men , 


a number of useful features have made it 


attractive fur developing expert systems. Followings are some of 
the important character isti cs of Turbo Prolog ' 

1. Turbo Prolog is descriptive 

Instead of a series of steps specifying how the 
computer must work to solve a problem, a Turbo Prolog 
programm consists of a description of the problem 
data. 

2 . Foreign language calls 

Every language has its own limitations. Turbo Prolog 
is also not an exception. FORTRAN is good for 
numerical computation, C and Assembly are handy for 
for lower level access of the machine. All this 
resources can be exploited by using Turbo Prolog's 
interface to C, Assembly, FORTRAN and PASCAL. 
Similarly Turbo Prolog routines can also be called 
from foreign languages. 

3. Graphics 

Turbo Prolog supports the full featured Borland 
Graphics, which has a strength of 60 different 
graphics predicates. 

4. Internal and external database 

Turbo Prolog's internal database can be used for 
storing informations in RAM at run time, which is 
limited by the hardware, usually upto 64GK bytes. 
External database is another feature which makes 
Turbo Prolog efficient for handling large amount of 
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data. Go ph i s t i c.ited schemes are available -Tor storing 
d a t a i ri a b i nat y t r e e . 


4.4 Parsing 

Parsing in an important part of liA. Since user is 
allowed to present the problem syirtbo J i cal ly a robust parser is 
needed to identify the problem. Once the string input is received 
a binary tree representing the equation is generated. A preorder 
search acotnpl i shes the equation identification job by pattern 
matching with pr e-dof ined templates. Terms are simplified by 
calling a predicate simp of arity 2. The simp clause is shown 
h ere : 

s imp ( In i t ial__l ist » Final„list) 

1 i st_exp( Initial__l ist , Initial__exp) , 
reduce ( Initial_exp,Fi nal_e:<p) , 
exp_l ist (Final_exp,Fi nal_l ist ) . 

Here predicate list_exp converts the string list to a recusive 
binary tree Initial„exp; reduce simplifies the binary tree 
Jnit i al_exp to Final_e;:p and exp_list converts a binary tree to a 
string list. 

Let us cusider an example: 

Expression to be parsed - (2-*uim1-4*ui+2-»uip1 ) /Dx 

At first this string is converted into a list of tokens. List_exp 
creates the binary parse tree of the list, which looks like- 

di v( plus (minus (mult( real (2> , var ("uiml" ) ) „ 
mult ( real (4) , var ( "ui") ) ,mult ( real (2) , var ( "uipl") > * var ( “Dx") ) 
Predicate reduce, which does the symbolic computation, will be 
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discussed in the following section. 


4.5 Symbolic Computation 

Symbolic computation is necessaiy as the form of the 
given equation is unknown and the functions are of unknown 
combination of the primitives. Symbolic computation is done in 
a recursive procedure, which handle strings in the form of binary 
tree. All tho mathematical functions sin, cos, tan, cosec, sec, 
cot. In, exp are available for use. MA has two symbolic 
computation modules one is written in Prolog and the other is in 
C. 1 h e main problem of calling procedures written in Prolog from 
C is that the heap memory used by these procedures is not 
released even after they succeed or fail. This is because, the 
C module controls the memory management while calling Prolog 
procedures. As a result the system soon runs out of memory wlti ch 
puts a severe limitation on rigorous number crunching. This 
problem is solved by developing a symbolic computation module in 
C, which handles a similar data structure, in case of nonlinear 
coefficients this module comes in much of use. 

We shall explore the reduce predicate, mentioned in the above 
section, now; reduce simlifies the expression by a preorder 
traversal of the parse tree. Look at the following primitive 
symbolic computation module, written in Prolog, which allows 
operators ", "/" only. 

/* Addition */ reduce ( plus ( X , Y) ,R> 

reduce C X , R 1 ) , 
r educe < Y ,R2) , 
plur <R1 ,R2,R> . 
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/* subs tract ion #/ 


/ si mu 3 1 i p I i cat j. on #/ 


/■si- J i v i. i> .i o n it' / 


/it default it / 


r e due e ( m i nus < X , Y > , R ) : 
reduce ( X , R1 ) , 
r educe < Y,R2> , 
mi nr (R1 , R2, R ) . 


r e du ce ( mul t < X , Y ) , R ) s - 
r educe ( X , R 1 ) , 
reduce (Y,R2> , 
mul r { R 1 , RH , R ) . 

reduce(div(X,Y) , R ) 
reduce ( X , R 1 ) , 
r educe ( Y , RS ) , 
d i v r < R 1 , R2 j R ) . 

reduce (R , R ) 


p 1 u r ( r ea 1 ( X ) ? r ea 1 C Y ) , R ) s — 

R--X+Y. 

m i n r ( r ea 1 ( X ) , r ea 1 ( Y ) , R ) :~ 

R-X--Y . 

mul r ( real C X ), real (Y) ,R ) !~ 

R=X-*Y. 

di vr ( real ( X ), real (Y) ,R ) :-~ 

R=X/Y. 

reduce calls itself recursively till its arguments become the 
primitive, » var < . ) or real(.>. Then plur, mirtr , mulr, divr do the 
algebraic, computation on the primitives. 

4.6 Handling boundary conditions 

MA treats the boundary conditions in a generic way. At 
first it identifies the type of boundary condition. In case of 
cau thy conditions the equation is discretized and simplified 
after substituting the grid parameters; discretized equation for 
internal grid points is equated with the discretized boundary 
condition for eleminating the fictitous grid point and the final 
boundary equation is sent to the C module through a recursive 
data structure. In case of elliptic problems special measures 
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should. be taken 'for the corner nodes. Corner nodes sharing 
neumann conditions arc equated to e'i iminate the fictitous grid 
p o i n 1 1 *> » 


4.7 Explanation 

One of the important features of expert systems is 
t ho capability of answering the user query. It has both 
psychological and material importance; it often give the user a 
feeling of assistance not as of being commanded by the system. 
Moreover system explanations are useful for a novice, who should 
be able to look through the system. 

The complete system is classified into several levels. A track of 
all the iule.» fired at every level is maintained. When the user 
typo why, system gets the recent rule and level from the database 
and answers it ’intelligently*. Although the technique for such 
quest i on- answer ing is not very sophisticated, for our purpose it 
does the job effectively, as the user session is comparatively 
simpler here. 
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Chapter Five - Conclusions 
5.1 Discussion 

As a very limited number of numerical export systems are 
available for persona] computer useis MA ran bo of much use m 
the scientific and engineering applications. MA can handle a wide 
range of RDEs with different kinds of boundary conditions. 
Presently available state- of-ar t analysis softwares CANSYS, 
NA5TFAN, IDEAS etc. 3 can handle only specific types of PDFs. MA 
certainly overcomes that weakness by addressing the problem in 
a generic way . The knowlwdge acquisition module of MA makes it 
dynamic, to let a numerical analyst experiment with new schemes. 
The unique Feature of handling directly the continuous physics 
problem is handy and approaches towards a natural way of 
thinking. The present prototype version of MA has several 
limitations too and those will be discussed in the following 
se ct ion. 


5.2 Limitations and Scope for future work 

As mentioned earlier, the major limitation of MA is that it can 
only deal with rectangular domain. Provisions are kept for 
interfacing MA with a solid rnodula which will expand the problem 
domain . Some of the features of symbolic computation module are 
ad hoc; for example, linearizing of the nonlinear source term may 
sometime generate even simplified expressions in a crude form. 
This is because it cannot put together complex similar terms. But 
this does not result errors in computation. So far the system is 
tested for a variety of problems. Graphics module needs more 
attention. At present MA is equiped with two dimensional 



gi'rtplri. 31) graphics is essential especially -for plotting 
sn ] u t j. ans of elliptic equations. Error handling file is 
incomplete; it needs a rigorous long system testing along with 
t, i i e u :> t» r f e e d La c K s . 

Conclusion 

A new problem solving environment, MA for solving PDEs is 
proposed and implemented. This work several unique features along 
with some drawbacks, discussed earlier. Earlier works on PDE 
solver were limited either by their too rigid handling of the 
problem (procedural subroutines) or by their inefficient number 
crunching power (systems written completely in relational 
languages). Although some expert systems for solving PDEs , like 
AGMEf! , F1DIS0L, Elliptic Expert , PLOD/PC C1-43 etc. are capable of 
tackling a wide range of PDEs, but all of them generate 
pro< oduraj codes, winch need to carry a procedural compiler 
(either internally or externally) and thereby demand more disk 
and memory space; moreover this is time consuming too. MA has 
overcome this difficulty. All oP the above mentioned systems 
address the problem as a PDE not as a phenomena. This is really 
an unwanted feature of an expert system, which is supposed to 
provide a higher level of user interface. MA certainly got rid of 
this limitation. Although the storage of continuous physics 
knowledge is somewhat ad hoc, this can be overcome by deriving 
governing PDEs from the fundamental equations using symbolic 
operators . 


48 











Fig.i 


>0 











| OBJECT 1 


r 


t 


EOUATION SPATIAL DOMAIN 

j r J— i r 


TIME DOMAIN 
I 



2D TIME DEPENDENT TIME DEPENDENT 


LINEAR NON-LINEAR 


I 1 1 

ELLIPTIC PARABOLIC HYPERBOLIC 


Fig. 3 


51 


CENTRAL 

1 i. T.. 


UE-RARY 

r '\;.PU 0 . 




- EXPERT SYSTEM 11ft- 
Craphics 



Quit 




Discretization scheae 

Conputdtioiul scheme 

OS shell 

Kelp 


f' 

* 



5Z 



FtG. 5 


GRAPHICAL OUTPUT OF MA 


53 



Fie. 6 FEEDBACK OF Mft T’RE. PROCESSOR 


84 


Welcome to EXPERT SYSTEM MA 
This is an expert system for solving 
Partial Differential Equations 

AU copyrights reserved to 
Hillol Kargupta 
IIT Kanpur, India 


Fig.7 - System plate no.1 
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Fig. 11 - System plate no. 4 
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Please, define the unknown term(s) 
found in the equation. Note that, 
unknown terms can either be a 
function of the dependent variable 
and independent variables or constants. 


diffusivity= 0.01 


• 1 

sc "Main menu] |? - Help] 


uuaaaa 


jwhy - Ex^laiTTl 


Fig. 12 - System plate no. 5 
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EQUATION TO BE SOLVED : 0.01*d2u/dx2-du/dt=0 

THIS IS A LINEAR, PARABOLIC IN SPACE, FIRST ORDER IN TIME, 
AND SECOND ORDER IN SPATIAL COORDINATES 

We need TWO boundary conditions along the X-axis and 
ONE initial condition to solve the problem. 

I will be asking you about the boundary St initial conditions 
and the grid parameters, soon. 

Press any key to proceed 


E sc - Main menu| |? - Help^ Quin jwhy - Explain j 


Fig. 13 - System plate no. 6 
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Enter X direction width =1 

Enter total number of X grid points, Nx =11 

Enter time step, Dt =0.001 

Enter total number of time steps, Nt =30 

Enter initial distribution, Uo =100 

First boundary condition =1000 

Second boundary condition =100 


IFsT^T^rnenul [?~-"HeTpl 


umirm 


Iwhy - Explain | 


Fig. 14 - System plate no. 7 












Fig. 15 - System plate no. % 
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6.1 Getting Started • 

This chapter describes the basic: operation of the expert system 
MA, including how to use the menu system and control the 
env i r onrnent . 

The README! file on the distribution disk gives a complete list of 
the files supplied on the distribution disks. MA comes ready 
install and run on an IBM PC or fully compatible computer. 

Follow one of these two procedures to load MA into your system s 

# If you are running MA from floppy disks ; 

1. Turn on and boot up your computer. 

2. Take your system start-up disk out of the disk 
drive. 

3. At the A; \>prompt , type ma and press Enter. 

# If you are running MA from your hard disk : 

1. Turn on your computer. 

2. Log onto the hard disk and change to the directory 
where MA resides. 

3. Type ma and press Enter. 

You should now see the logon message shown in Fig. 7. Press any 
character key; MA’s main window and menu will appear, as shown in 
F i. g . 8 . 

6.2 A Quick Guide to the Menus and Hot Keys 

In this section an introduction to MA's menus and some special 
shortcut keys (known in the vernacular as hot keys) is given. 
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The Main Menu : 


MA ' y main menu showy you the commands and pulldown menus 
available. You can select an item from the main menu in one of 
the Following two ways s 

1. When main menu is active, by pressing the associated 
highlighted capital letter (S for Solve, 0 for Options, 
etc . ) 

2. When main menu is active, by first moving the highlight 
bar with the arrow keys, and then pressing Enter 


The Pull-down Menus: 

Four of the items on the mam menu invoke pull-down menus when 
you press them; these pull-down menu items are Solve, Utilities, 
Database and Options. Here’s a brief description of the main menu 
items : 

a) Solve 


Using this item you can invoke the problem solving 
module, which initiates a course oF dialog with the user 
and finally solves the problem. Two options are offered, 
namely s 

Feed equation - You can directly feed the PDE to 
be solved using this facility. 

State problem - The physical problem can Instated , 
which in turn fetches the PDE and relevant informations 
from the knowledge base. 



b> 


Utilities 


The items on this menu allow you to handle •Tiles, 
manipulate directory , invoke DOS, call graphics module. 

I he listing oT the items are s 

1. Copy file One can copy a file to a destination file by 

using this option. 

?. Print file This can be invoked to get the printout of a 

file directly from the environment. 

3. OS shell One cati go to the operating system directly 

from the system; type exit to come back 
to the system. 

4. Graph Plot One can directly call the graphics module 

for plotting graphs, by stating the 
graphics filename. 

5. Sketch Pad At present no interface is provided 

with MA. This is kept for future 

modifications, which will enable the 
user to define the problem domain 

graph i cal ly . 

c) Database 

Database option handles the manipulation of the 
knowledgebase. Four choices are available: 

1. Display database file - All the updatable knowledge 
is stored in 1 1 » e file mabase.dba. This option allows the 
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u'*«m ^ o look at the da t abas e without allowing any 
change. Press Esc key for coming out from display mode. 

2. Update database -■ System knowledgebase can bo 

updated using this option. At first system will ask you 
about the physical name of the problem, eg. heat 
conduction, vibration etc. and then the governing PDE . 
Variable properties can be introduced. Use default 
denotations for variables, otherwise define new 
variables using menu Opt ions/Var iable denotations. 

Remember once you update the database all the new 
informations are kept in RAM only, not in the file. So 
you must invoke Save database to file menu if you want 
to retain it. 

3. Edit database file - Any bug in the database can be 
eleminated using this option. This option invoke an 
editor loaded with the database file "mabase . dba" . 

4. Save database in file - This option saves the 
current database in the file "mbase . dba" . It always 
creates a back-up, named mabase.bak, 

d) Options 

This menu provides several subsequent pull-down menus, 
from which you can change the default description of the 
dependent and inde pen tent variables, computational 
scheme and discret ization scheme. 



"I. Variable denotation 


Different denotations for different 


type of problems makes the system handy 
and flexible. Default denotation for the 
dependent variable is u, x-axis x, and 
y-axis y. 

H. Discret izat ion scheme User can specify the discretization 

scheme beforehand. This can also be dune 
during the dialogue session with the 
system. At present Discretization scheme 
offers a choice among the Fully implicit 
scheme, Crank-Ni col son scheme and 
explicit scheme. 

3. Computational scheme Similarly the computational scheme 

can also be beforehand. At present 
Succesive Over Relaxation (SOR), Gauss- 
Siedel, Tri-diagonal algorithm, and 
Gauss-Jordan methods sire available for 
solving the set of algebraic equations. 


e> Quit 

This item should be chosen to come out of 

the system. 


Hot Keys and commands 

esc. go to main menu 

? invoke on line help file 

why query about the current process 

quit 

?'•> 
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6.3 Solving a Partial Differential Equation s 


Solve option provides two choices for solving the 
problem. The State problem option invokes the 
preprocessor directly and displays a pop-up menu, 
containing all the physical problems, defined in the 
database, and then asks the user about the unknowns 
present in the equation, if any. Fig. 9 shows the 
procedural flow of HA. 

If the Feed equation option of the main menu is hit, system 
will want the equation to be solved, by prompting "Enter 
the equation to be solved", as shown in Fig 10; after 
parsing the equation the list of unidentified terms are 
sent back to the user for proper defination. Once the 
equation is completely defined the processed equation is 
sent back to the user for approval along with several 
useful informations about it. This is shown in Fig. 6. 
If the discretization scheme is not chosen beforehand a 
menu offering several discretization schemes appears in 
the screen. After this the dialog session for obtaining 
the grid parameters and boundary conditions comes into 
action. After the boundary conditions are identified and 
processed, discretized coefficients are sent ba c k to the 
user. Then system suggests the best computational 
scheme; but the user can deny' it and choose the desired 
one from the menu, as shown in Fig 6. At any stage 
user may query the system by typing why. Now the control 
goes to the procedural number crunching modules. Output 
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re<iu3t is stored m both an output file and a graphics 
data -File. Results are also displayed on the screen at 
the same time. User can stop the number crunching 
process by pressing Ctrl— break key. As soon as the 
computation work is over message window prompt will let 
it know. A pop-up menu appears in the mid-screen which 
offers several actions. Options are equation dependent. 
Followings are the different options available for 
different type of equations s 

a. Elliptic equation 

1) U (dependent variable) vs. X (x-axis variable) 
If this option is chosen the y-values for which a U/X 
graph will be plotted, are further asked. 

2) U (dependent variable) vs. Y (y-axis variable) 
If this option is chosen the x-values for which a U/Y 
graph will be plotted, are further asked. 

3) Dad- to main menu 

b. Parabolic equation 

1) Simple plot - This option offers an usual graph 
plotting with the spatial domain along x axis and the 
dependent vaiiblo along y-axis. 

2) Back to main menu 
b . Hy pe r bo 1 i c e q ua t i on 

1) Simple plot - This option offers an usual graph 
plotting with the spatial domain along x axis and the 
dependent varible along y-axis. 


2) Back to main menu 



6.5 Mathematical functions 


Several mathematical functions are available in 
the library. A listing of these functions are 5 


s i n ( X > 

sm function 

c.os(X> 

cos function 

tan ( X ) 

tan function 

cosec ( X ) 

cosec function 

se c ( X ) 

sec function 

cot < X ) 

cot function 

In ( X ) 

natural logarithm of X(base e) 

1 o g < X ) 

base 10 logarithm of X 

e x p ( X ) 

e raised to the value of X 

abs(X) 

absolute value of X 

X"Y 

X to the power Y 


6.6 Errors 

A number of error* may occur at run time. All the 

errors along with their reasons are stored m MA.ERR. In case of 
any error system sents a message to the message window. If the 

user commits any mistake in typing either the aquation or the 

paramo tors , system points out the mistake immediate ly and lets the 
user correct it. In case of memory overflows, which may occur due 
to extensive grid resolution, very large number of iterations etc., 
system sends a message in the message window and resets- the 

system. 





Appendix B — Case study 


Expert system MA has a robust FDM solver, which can handle a 
wide range of partial differential equations. Our system has 
been tested for a wide range of PDFs and so far at has proved 
its computational ability. In this section we shall discuss a 
sample problem, to illustrate the complete process of solving a 
PDF using MA. Although a simple problem is taken here, MA can 
handle complex situations too. 

Consider the following heat transfer problem : 

a) The wall of a furnace is 1ft thick and is made of brick, which 
has a thermal diffusivity of 0.01 ft a /h. The temperature of the 
wall is 100 °F when the furnace is off. When the furnace is 
fired, the temperature on the inside face of the wall reaches 
1000°r quite rapidly. The temperature of the outside face of the 
wall is maintained at 100°F by natural convection. Determine the 
evolution of temperature profiles within the brick wall. 

b) Insulation is placed on the outside surface of the wall. 

Assume this is perfect insulation and show the evolution of the 
temperature profiles within the wall when the furnace is fired to 
•1000°F. 

Both the above problems can be solved by MA 
with ease. Fig. 11 shows the problem type menu, which comes as 
soon as one hits the State problem option. Heat conduction (ID) 
can be chosen, which results in the preprocessing session, as 
shown in Fig. 12, The System then asks for the value of diffusivity. 
Once all the terms of the PDF, governing the problem are defined 
a feedback message is sent to the user. This is shown in Fig. 13. 
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Now the processor module is invoked, which asks the user about 
the grid parameters and boundary conditions as shown in Fig. 14. 
All the values are entered as stated in the -first problem; 
finally the computation module is invoked, which saves the output 
to the user specified file, displays on the screen and also keeps 
in a database chain for graphics. Graphics output of the first 
problem is shown in Fig. 16. 

Fig. 15, Fig. 17 show the changed boundary condition and the 
co r res ponding solution in graphics, respectively. 



